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Nonlocality can be quantified by the violation of a Bell inequality. Since this violation may be 
amplified by local operations an alternative measure has been proposed - distillable nonlocality. 
The alternative measure is difficult to calculate exactly due to the double exponential growth of 
the parameter space. In this article we give a way to bound the distillable nonlocality of a resource 
by the solutions to a related optimization problem. Our upper bounds are exponentially easier to 
compute than the exact value and are shown to be meaningful in general and tight in some cases. 
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I. INTRODUCTION 

When two separated parts of a quantum state are mea- 
sured in different bases, then the outcomes can be cor- 
related in a way that cannot be explained by informa- 
tion shared before the separation This property has 
been termed quantum nonlocality. Nonlocal correlations 
have since proved to be a useful resource leading to new 
applications such as, for example, quantum key distri- 
bution Stronger correlations that are still in ac- 
cordance with the nonsignaling postulate of relativity @ 
can be defined and are formalized and studied in so-called 
generalized nonsignaling theories [3, Q in which quantum 
correlations are a special case. Assuming nonlocality that 
is super-quantum to some extent has interesting conse- 
quences for nonlocal c omp utation [9J and for communi- 
cation complexity [ToL lllj . Maximal nonlocality allows 
to compute every distributed Boolean function with just 
one communicated bit [l2j . 

Local correlations obey certain linear constraints, so- 
called Bell inequalities [1|. The extent by which a Bell 
inequality is violated by a correlation can be taken as 
a measure for nonlocality (see [Tijl - fTEj for other nonlo- 
cality measures and related results). The question arose 
whether this measure is also meaningful for quantifying 
the nonlocality of a resource consisting of several copies 
of a correlation. In a context where nonlocality is more 
useful the stronger it is this measure is problematic if 
stronger nonlocality can be obtained from a number of 
weakly nonlocal correlations. In other words, we want to 
know if two parties having access to correlations violating 
a Bell inequality by some small extent can execute local 
operations to obtain a higher violation. The two parties 
can carry out arbitrary local operations, but cannot com- 
municate. This process is called nonlocality distillation. 

It has recently been shown that the Clauser-Horner- 
Shimony-Holt inequality [16| (CHSH) as a measure for 
nonlocality in minimal dimensions is indeed problem- 
atic because it is distillable in general fl7| . A proto- 
col by Brunner and Skrzypczyk manages to dis- 
till an arbitrarily weak nonlocal correlation, that is still 
super-quantum, to the extent where communication com- 
plexity collapses. As proposed in 17] a more mean- 
ingful measure, which is by definition undistillable, is 
the maximal CHSH violation achievable from many re- 



alizations of a given nonlocal resource by any distilla- 
tion protocol. Brunner et al. [IH recently compared 
this alternative measure, termed distillable nonlocality, 
to the Elitzur-Popescu-Rohrlich (EPR2) decomposition 
approach 0, HH ■ The authors discovered examples of 
bound nonlocality and activation - the box-world ana- 
logues to bound entanglement and its activation. 

H0yer and Rashid [l9| thoroughly analyzed differ- 
ent distillation protocols, proved optimality in restricted 
classes and discovered a strong dependence between the 
optimal protocol and the parameters of the resource that 
it distills. Today our knowledge about general distilla- 
tion protocols is still very limited. Although tools and 
techniques to analyze fixed protocols exist (for example 
by discrete maps [IHHil), difficulties arise and the lack 
of a neat framework becomes obvious when attempting 
to calculate the distillable nonlocality of a fixed resource. 
Besides what results from some trivial symmetries not 
much is known to simplify the task. It usually requires an 
exhaustive search over all possible distillation protocols. 
The fact that the search space grows doubly exponential 
in the number of involved correlation copies represents 
a serious problem for the usefulness and understanding 
of this alternative nonlocality measure. Therefore, as re- 
cently suggested in [l8[ , meaningful bounds for distillable 
nonlocality that can be derived more efficiently are a rea- 
sonable compromise. 

In this article we propose such bounds and analyze 
them. They rely on a recursive optimization problem 
that has a strong connection to the process of optimiz- 
ing distillation protocols for isotropic nonlocality. Solved 
with a dynamic-programming approach our bounds yield 
an exponential gain in the run-time compared to a brute- 
force search. The efficiency advantage enables one to 
analyze instances with the number of available correla- 
tions limited by 9. To demonstrate this we calculated the 
distillable nonlocality bounds for a set of fixed isotropic 
systems and found that, in this context, our bounds 
rule out distillation and are therefore tight. Further- 
more, we present a general idea how to extend isotropic 
bounds to the non-isotropic case. Here, our calculated 
bounds show that the distillable nonlocality approaches 
the CHSH nonlocality of a single copy when the resource 
is chosen closer to an isotropic line. 
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II. PRELIMINARIES 

Following a general approach to nonlocality |8j , we con- 
sider correlations in the joint behavior of the two ends of 
a bipartite input-output system, characterized by joint 
probability distributions PVL on random variables A, B 
for each x,y. Let x and a be the input and output on 
the left-hand side of the system, and y and b the corre- 
sponding values on the right-hand side. On inputs x and 
y the system returns outputs a and b with probability 
P x AB {a,b). 
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In a distillation setting two parties, Alice and Bob, 
sharing a number of systems can independently choose 
inputs - possibly derived from local outputs of other sys- 
tems - and collect outputs on their ends of the systems. 
In other words, they can apply any classical circuitry 
to their local parts of the systems. Such a local input- 
output strategy is called a wiring l2l| . A party re- 
ceives its output from a system immediately after giving 
its input, independently of whether the other has given 
its input already. This prevents the parties from signal- 
ing by delaying their inputs. On the other hand it allows 
for wirings that have a different temporal order on the 
two sides. For example: While Alice's input into sys- 
tem 2 depends on the output of system 1, Bob's input 
into system 1 depends on the output of system 2. If the 
temporal order on both sides is equal we call the wiring 
ordered and disordered otherwise. In this paper we con- 
sider both kinds, summarized under general wirings as 
defined in [22| . Since we are only interested in optimal- 
ly ignoring non-deterministic strategies is sufficient due 
to convexity. Thus, we are allowed to ignore randomness 
shared between the two players to simplify the proofs. 

Here, we assume nonlocality in its simplest form, 
namely, in the binary setting where each input and each 
output has two possible values, i.e., a, b G {0,1} and 
x,y e {0,1}. We refer to [|[ for a detailed description 
of the convex polytope of binary nonsignaling systems 
(AfS). In the case where both inputs and both outputs 
are binary, the only Bell inequality (up to symmetries) 
is the CHSH inequality. Furthermore, the set of eight 
CHSH inequalities is complete for binary systems in the 
sense that if none of them is violated, then the system is 
local. A system is called isotropic if it is a probabilistic 
mixture of opposite nonlocal vertices of AfS, which are 
known as PR-boxes [f|. See [23j for a formal definition 
of isotropic systems. 

Let P abbreviate any binary system. Suppose that 
Alice and Bob share n copies of P, which we index as 
Pi,...,P n . The outputs of the i-th system Pi are de- 
noted (ai,bi) <E {0, l} 2 . So, the binary strings a = 
(oi, 02, a n ) G {0, 1}" and b = (b u b 2 , &„) G {0, 1}™ 



are the outputs of the n systems to Alice and to Bob. 
Also, from now on, let A = (Ai,...,A n ) and B = 
(Bi, B n ) denote random variables for the strings of 
the n collected bits on both sides. Any wiring of n shared 
copies of a system then fully determines a joint distribu- 
tion Wab on the space of all sequences (a, b) G {0, 1}™ X ". 
Alice and Bob can condition their wirings on the inputs 
x and y. We denote a system of wiring distributions on 
n copies of the binary system P by P n . We will use the 
matrix representation 



pn _ 



W AB 



where WVg, for any x,y G {0,1}, is a 2™ x 2™ matrix 
with the probability W^ v B (a, b) at position (a, b). As the 
final step in a distillation attempt the two parties map 
their inputs x and y and the collected strings a and b to 
one local bit each. For all x, y let f x ,g y ■ {0, 1}™ — > {0, 1} 
stand for these local Boolean functions. It will often be 
convenient to write them as truth tables fxigy G {0, l} 2 . 
Now we have all the ingredients to define a distillation 
protocol on binary systems. 

Definition 1. A deterministic nonlocality distillation 
protocol on n copies of a system P , denoted DP = 
(P n , /o, fx, go, 9i), consists of a system of wiring distri- 
butions P" and local decision functions /o,/i and go,gi 
for both inputs and both sides. 

We refer to [l6| for an exact formulation of the CHSH 
nonlocality and to [l7] for the definition of the CHSH 
nonlocality of a single binary system NL(P) as required 
here. Note that NL(P) > 2 indicates that P can be 
used to violate a CHSH inequality and is therefore called 
nonlocal. To measure the violation of a CHSH inequal- 
ity by a distillation protocol we use a slightly different 
representation. Observe that for each input pair x, y the 
correlation function of the binary system simulated by 
a distillation protocol can be described by the following 
inner product: 

(f x ,g y ) = E[f x (a) = g y (b)} - E[f x (a) + g y (b)\ 

= X)W^(a,6)(l-2/ a (a))(l-2 5l ,(6)) 

ab 

= (l-2/«(l-2s,). 

The two expectation values are calculated over all pairs 
(a, b) ~ W^b- The nonlocality of a protocol is the max- 
imal violation of a CHSH inequality by the simulated 
system. 

Definition 2. The CHSH value of a distillation protocol 
DP = (P n , fa, fi,go,gi) is measured by 

NL(DP) = max|(/„ 5y ) + (f x ,g y ) + (f x ,g y ) - (f x ,g y )\ 

xy 

where we use x and y to indicate bit flips, that is, 
= 1 and 1 = 0. Note that NL(DP) > 2 indicates 
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that the protocol simulates a nonlocal system, whereas 
NL(DP) > NL{P) indicates a successful distillation at- 
tempt. 

We will henceforth assume that the maximum of this 
expression is reached at xy — 00. Our reasoning can 
easily be extended to all four cases. 

Motivated by the discovery of successful distillation 
protocols for CHSH nonlocality another measure has 
been proposed that expresses the nonlocality of the op- 
timal protocol on a given resource. Brunner et al. [18[ 
have recently formally defined and analyzed distillable 
nonlocality in comparison with the EPR2 decomposi- 
tion [Til, HU ■ Here, we use a slightly different notation 
that expresses the same idea. 

Definition 3. The distillable nonlocality of n copies of 
a system P is defined as 



D(n,P) 



maxNL(DP), 

DP 



(1) 



where we maximize over all deterministic nonlocality dis- 
tillation protocols on n copies of P. 

For determining distillable nonlocality precisely by 
an exhaustive search over the space of all possible se- 
quences (P n , fo, /i, go, gi) one requires to test roughly 
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instances, calcu- 



lated as the product of all involved Boolean functions. 
Because this is infeasible with normal hardware already 
for n > 2, giving an asymptotic bound to D(n, P) instead 
which can be computed more efficiently, is a reasonable 
compromise. 



III. RESULTS 

First, we show an upper bound on D(n, P) if P is an 
isotropic system (Corollary Q]) and move on to a general 
upper bound afterward (Corollary [5]). We start by ma- 
nipulating (JXJ) to reveal the core of the task. One can 
group the Boolean functions fo, fx, go, gi, over which we 
optimize in (JlJ , by their output distribution. For a given 
number n, a class Ck of Boolean functions is defined as 

C fc = {/:{0,ir-»-{0 J l}:|/- 1 (l)| = fc}. 
We can now rewrite the distillable nonlocality as 
D(n,P)= max maxNL(DP), 

0<fc x ,;a<2" DP 

where for each k x ,l y we maximize over all distillation 
protocols restricted by f x £ Ck* and g y £ Civ for all x, y. 
Now concentrate on maxDpNL(DP) with fixed k x and 
l y for all x, y. According to Definition [2j this subproblem 
boils down to finding the system of wiring distributions 
P n and functions /o £ C fc o,/i £ C k i,go £ C ; o,5i £ C; 1 , 
such that the term 

(Jo, so) + (A, go) + (/o,.gi) - (h,gi) 



is minimal/maximal. Since for any x,y we have 
(f x ,g v ) = (l-2f x ) T Wl y B (l-2g y ) 



1 - k x /2> 



iy/2 n -'+4fJW* A %g y , 



bounds for the correlation functions (f x ,g y ) can be de- 
rived from optimizing f T WABg, with Was a wiring dis- 
tribution and / £ Cfcx , g £ Civ , independently for all 
x,y. We will derive bounds to / t W / ab5 by slightly re- 
laxing the constraint that Wab needs to be a distribution 
obtained from a wiring. This relaxation has the following 
background: 

Suppose fixed functions / £ C'k* , g £ Civ and a fixed 
wiring distribution Wab- Let (ai,b±) be the outputs of 
the first system Alice accesses in this wiring. The k x 
preimages of 1 under / are split into two parts with ref- 
erence to the output bit a\ and the l y preimages of 1 
under g are split into two parts according to b\, i.e., we 
have 



and 



fc 3 '=|{/- 1 (l):a 1 =0}| + |{r 1 (l):a 1 = l}| 



iy=\{g-\l):b 1 = 0}\ + \{g- 1 (l):h=l}\. 



The wiring Wab provides a probability distribution on 
(ai,&i), possibly conditioned on outputs of some other 
shared systems. Roughly speaking we will find the opti- 
mal splitting of /, g and the optimal distribution for the 
output pair (oi, &i) by considering the optimal splittings 
and optimal distributions in the four subproblems, where 
the pair (ai,&i) is assumed to have fixed values. In this 
way we obtain a recursive n-level hierarchy of optimiza- 
tion problems. The following definition makes this idea 
precise. 

A problem instance is given by the tuple (P,n,k,l), 
with a binary system P, integers n > 1 and < k, I < 2". 
The recursive optimization to solve is: 



Definition 4. Given any instance (P,n,k,l) fix the pa- 
rameter p = Pab(®i 0) and find the values 



S+(k,l) = 



S-(k,l) 



Hh-pM-i(i,l-3) + S+-i(*-i,3)], 
ifn>0, 
kL otherwise, 



iny p[S n _ 1 (i,j) + 6 n _ 1 (k - i,l - j)] 

+(2-p)[*n-i(«»'-i) + *n-i(*-*.i)] 5 

ifn>0, 
kl, otherwise, 



where we optimize over the ranges k — min(fc, 2™ ) < i < 
mm(k,2 n - r ) and I - min(Z, 2 n ~ x ) <j < min(Z, 2™" 1 ). 
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Solutions 5+(k x ,l y ) and 5~(k x ,l v ) to the instance 
(P, n, k x , l y ) capture the possibilities we have when con- 
structing a wiring distribution Wab and functions / e 
Cfcx , g £ Civ , such that the product f T WABg is optimal. 
Since in the above optimization the subproblems on the 
same level are handled independently we get better so- 
lutions than with distributions from real wirings, where 
this independence is not given. 

Lemma 1. Let P be any isotropic system and suppose 
integers < k, I < 2™. For any wiring distribution Wab 
on n copies of P and functions f £ Cfc, g € Ci we have 

S-(k,l)<f T WAB9<S+(k,l), 

where S~(k, I) and <5+(fc, I) are the solutions to (P, n, k, I). 

Proof. The idea for proving the upper bound is to decom- 
pose Wab into four subwirings and /, g into two func- 
tions each, on which we then inductively prove the state- 
ment. The decomposition is done with reference to the 
first system Alice accesses, called Pi, its outputs denoted 
(oi, b\). We write a\ for the string (<22, a n ) and b\ for 
(62, b n ) and accordingly we use the random variables 
Aj = (A 2 , A n ) and B 1 = (B 2 , B n ). 

We fix the probability p = P AB (0, 0). Let the Boolean 
function h : {0, l}™ -1 — >• {0,1} determine a bit depend- 
ing on Bobs local outputs b\ — (62, ■■■,b n ) of the n — 1 
systems P2, ...,P„ shared with Alice, such that 

W AlBl]Bl (0,h(b 1 )\b 1 )=p 

for all b\ £ {0, l} n_1 . Note that if the wiring defining 
Wab is ordered then Alice's and Bobs inputs into the first 
system are constant and thus h is also constant. Other- 
wise, only Alice's input is guaranteed to be constant. In 
this case h is used to react to possible changes caused by 
Bobs input into the first system, h always exists since 
isotropic systems have symmetric output distributions, 
and therefore for any 6j the first system Pi outputs ei- 
ther (0,0) or (0, 1) with probability p. 

According to the result of h we split Wab into four 
parts by grouping all probabilities Wab (a, b) by four 
possible outcome pairs of the first system. In the same 
manner we split the functions /, g into the subfunctions 
f'( a l) = /(0,ai) and f"(ai) = /(Mi) for any aj and 
9'{h) = 9{Kh), h) and g"{bi) = g{h{b%), h) for any b t . 
Applied to f'W ABg we can identify four parts as 

f T W AB g = f( a ) W AB(a, b)g{b) 

ab 

= [f{ai)WAB{QaiMh)bl)g'{h) 

ajbj 

+ f"(a l )W A B(la l Mh)bi)g"(b l ) 
+ f(a l )W A B(la l7 h(b l )b l )g'(b l ) 

+ f(a- 1 )W A B(0a- 1 Mh)h)g"(bi)]. 

Now observe the probability Wab(0oi, h(bi)bi). By the 
definition of the function h we can derive 

WA^&hih)) = J2WB 1 (bi)W Al B 1 \BMKk)\k) =P 

h 



and therefore, if we factor out the distribution of the first 
output pair (ai, 61) we get 

W A B(0ai,h{bi)bi) = pW A - 1 b- 1 \a 1 b 1 &l|0, h(bi)). 

For any oj, b\ the probability W Ai b i \AiB 1 (oi, &l|0, h(bi)) 
is independent of the distribution of the first output pair. 
That means Pi can as well be replaced by local circuitry. 
An alternative wiring distribution is obtained from the 
original wiring by fixing ai — and &i = h{b\) for 
any aj,&j. Therefore, there exists a wiring distribution 
Wa 1 b- 1 on n — 1 systems, for example the one described 
above, such that for each a%,bf we have 

W AB {Qa 1 ,h{h)h)=pWA l B l {a-M- 

The above reasoning can be applied to all four parts of 
f T WABg implying the existence of wiring distributions 
W A % , W A % , W$ Bl and W A % , such that 

fw AB g = p (r<v + r T <V) 

+ (l/2-p) (f" T W A % l9 ' + f' T W^ Bl g") . 

Assume now the induction hypothesis that for any wiring 
distribution Wa- x b- x on n — 1 isotropic systems and func- 
tions /, g : {0, l} n_1 — > {0, 1} we have the upper bound 

fWA-^g^S+^Qf-HVllg-^l)]). 

Since the basis f T P AB g < S+G/ -1 ^)!* Iff"^ 1 )!) trivially 
holds, we can conclude 

fW AB g<pS+_ 1 (\f'- 1 (l)\,\g'- 1 (l)\) 

+ pS+-i(\f"- 1 (i)\,\g"- 1 (i)\) 

+ (l/2-p)Ci(l/'~ 1 (l)l.lfl"" 1 (l)l) 

+ (i/2- P )t 1 (ir i (i)i ! r i (i)i) 

We used that < < min(fc,2"- 1 ) 

and < Iflf'-^l)!,^"- 1 ^)! < mm(l,2 n - 1 ) holds. With 
the same arguments one can prove the lower bound 
f T W A sg > 5~ (fc, I) by induction. □ 

The shown bounds to f T W AB g yield the following 
bound to the distillable nonlocality of isotropic systems: 

Corollary 1. Let Pi SO be an isotropic system; then 

D(n, P iso ) < 2 - + i[SUk°, 1°) 

+ d+(k ,l 1 )+S+(k\l°)-d-(k\l 1 )}. 

where 5^(k x , l v ) and 8~{k x , l y ) are solutions to the prob- 
lem instance (Pi so ,n,k x ,l y ) for all x,y. 

Proof. Note that from the definition of the inner product 
it follows that it is sufficient to consider positive values 
to find the maximum (see Lemma [2] in Appendix (XJ . 
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Therefore we can omit the modulus in the CHSH value 
and get 

D(n,P) = max max NL(DP) 

0<fc* 2«<2™ DP 



max max 

0<k",ly<2 n DP 



(fo,9o) + (h,9o) 
+ (fo,9l) - (fi,9i) 



where the second maximum is over f x E and g y E Civ 
for all x,y. To complete the proof use (f x ,9y) = 1 — 
k x^ 2 n-i _ p/ 2 «-i + 4:f^Wl v B g y for all x,y and apply 
Lemma [TJ □ 

From the bound on distillable nonlocality of isotropic 
resources we derive a general bound as follows. 

Corollary 2. For any P, we have D(n,P) < D(n,Pi SO ) 
where Pi SO shall be isotropic and minimize NL(Pi so ) un- 
der the constraint P = qPi so +Pl(1 — q) for any q E [0,1] 
and any local system Pl ■ 

Proof. By contradiction. Assume n P's can be distilled 
to a value higher than D(n,Pi SO ). Given n P^o's we 
can construct a distillation protocol as follows: Combine 
each of the n Pj so 's with one local system Pl to get the n 
mixtures P — qPi SO + (1 — <z)Pl- Since we assumed that 
n P's can be distilled above D(n, Pi SO ) we constructed a 
protocol on n P^o's that achieves the same. □ 

Bounding D(n,Pi SO ) by Corollary [T] yields an expo- 
nential speedup compared to the naive search through 
all protocols. The idea is to use dynamic programming 
to solve the problem in Definition |4] in run-time 2°^ 
(Lemma[3]m Appendix [Bj . Bounding D(n,P) by Corol- 
lary [2] requires finding the least nonlocal isotropic sys- 
tem Piso, such that P can be obtained from a convex 
combination of Pi SO and a local system Pl. This can be 
done efficiently using a linear program (see Lemma [4] Ap- 
pendix [C]). Finally we calculate a bound on D(n,Pi SO ) 
as described above to complete the bound to D(n, P). 

To compare the results, an exhaustive search to de- 
termine D(n, P) exactly requires an effort of 2 2 ° {n) - an 
exponential increase to our solution. 



IV. APPLICATION 

Here, we illustrate Corollaries [T] and [5] on example ap- 
plications. For this purpose consider a two-dimensional 
section in the binary nonsignaling polytopc defined as 
the convex combination of a nonlocal vertex Pnl (a PR- 
box), an unbiased mixture of two local vertices Pc (per- 
fectly correlated random bits) and the local isotropic 
system Pp (l ying on the CHSH facet corresponding to 
Pnl) [II EMM ■ We mix as: 

P s ,5 = eP NL +6P C + (!-£- S)P F , 

where s S [0,1], S E [0,1] and e + 5 < 1 to ensure 
nonsignaling. By setting <5 = we obtain a set of isotropic 



systems P lso (e) = P e>0 , where NL(P iso (e)) = 2(e + 1). 
These are special when it comes to distillation because 
any system can be turned into an isotropic system while 
preserving nonlocality (a procedure known as depolariza- 
tion [13]). It is a long standing conjecture that isotropic 
systems cannot be distilled. Strong evidence to support 
this has been provided by Dukaric and Wolf [23], who 
showed that in the quantum region infinitely many ex- 
amples of asymptotically undistillable isotropic systems 
must exist. Also, Short [25[ gave a proof for two-copy 
impossibility, i.e., he showed that D(2,Pi SO ) = NL(Pi SO ) 
holds for any isotropic system. With the new bound 
at hand we can deliver additional evidence. We found 
that our isotropic bound (Corollary [T]), calculated ex- 
actly for rational instances with Maple [26], confirms 
D(n*,P lso (e)) = 2(e+l) = NL(P lso {e)) if walking on the 
isotropic line with step length a (e E {a, 2a, 1 — a}) 
in all tested cases: 



a 


10~ 6 


10~ 5 


10~ 4 


10~ 3 


10~ 2 


10" 1 4- 1 


n* 


3 


4 


5 


6 


7 


8 9 



An immediate consequence is that in these settings our 
bound is tight, the optimal protocol simply imitates 
Pi SO {e). Interestingly, the bound on D{n 1 P iso {e)) drops 
significantly below NL(Pi SO (e)) if we restrict the pro- 
tocol to use unbalanced functions, i.e., if we demand 
l/o^WU/f^ 1 )! 2™" 1 for example. Figure [I] illus- 
trates the isotropic bound and this effect by an example 
calculation. 




k° + k 







p + l 1 



FIG. 1. Plotted are the calculated bounds to D(6, -Pi/5,0) ( z ~ 
axis) for < kP+k 1 < 2 7 (as-axis) and < P+l 1 < 2 7 (t/-axis). 
The peak in the center is the CHSH value of a fully balanced 
protocol, i.e., = Iff^M] = 2 5 for all x,y. It's the 

only class that reaches 12/5 - exactly the value A r L(P 1/ / 5 o). 
Trivial protocols that ignore 5 of 6 systems achieve this value. 

Departure from the isotropic line yields another obser- 
vation based on the calculated isotropic bounds. Let e 
be fixed. By combining P £) o in a convex way with P £) i_ e 
- a system of the same nonlocality as P e .o sitting on the 
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nonsignaling facet (also called correlated NLB [HI, [III) ~~ 
we get the system 

P q = (l-q)P E , +qPeA-e- 

The bound in Corollary [5] allows more distillation on P q 
the closer the system approaches the nonsignaling facet. 
Starting with q = and D(n, P E<0 ) - NL(P efi ) = the 
difference between the bound and NL(P q ) grows with 
increasing weight q until D(n,P e ,i_ e ) < 4 - an asymp- 
totically tight bound for any P £: \- £ with e > as con- 
firmed by the distillation protocol in [ll|. Therefore, our 
tests support a suspected property of distillation pro- 
tocols, namely that, geometrically speaking, distilling a 
system pushes it towards the isotropic line. The longer 
this distance is, the better the system can potentially be 
distilled (See Fig. [5]). This is a property inherent to all 
bipartite nonlocality distillation protocols [ll], Qjl EH Hfll 
known to the author. 



bound for: 




p CHSH inequality p ( 



FIG. 2. Here, we illustrate the idea behind the general distil- 
lation bound of Corollary [2] in a two-dimensional wedge. The 
distillable nonlocality of P q (for q — 1/5,2/5,3/5,4/5,1) is 
shown to be limited by the bound to the weakest isotropic 
system Pi SO , such that P q can be expressed as a convex com- 
bination of Pi SO and Pc- 

Another implication of the calculated isotropic bounds 
and Corollary [5] is the existence of many sets closed un- 
der (n < 10)-copy wirings. We have established the fol- 
lowing: Given the set S(Pi SO ) of systems that are con- 
vex combinations of a certain Pi SO and any local system. 
Then D(n,Pi SO ) — NL(Pi SO ) implies the existence of a 
set S'(Pi SO ), where S(Pi SO ) C S"(P iso ), that is closed un- 
der n-copy wirings of systems in S(Pi SO ). 

Note that the consequences of Corollary [5] are in- 
dependent of how the isotropic bounds are derived. 
Since the corollary expresses a generic dependency be- 
tween isotropic and general distillable nonlocality any 
isotropic statement has direct implications in the gen- 
eral case. Therefore, we can use known asymptotic 



distillation bounds to generalize the above reasoning: 
Since D(oo,P iso ) = NL(P iso ) holds for infinitely many 
isotropic systems in the quantum region 23], the exis- 
tence of infinitely many sets, each with a different nonlo- 
cality limit, that are closed under wirings follows directly. 
This answers an open question posed in [20| . 

V. CONCLUSION 

Distillable nonlocality as a measure for the nonlocality 
of a given set of correlations is sometimes preferable to 
the original CHSH nonlocality, which can be amplified 
by local operations and is therefore problematic in cer- 
tain settings. Here, we derived bounds for the alternative 
measure. Bounding distillable nonlocality is useful if the 
bound can be calculated more efficiently than the exact 
value and is tight enough to be meaningful. We show 
both properties: An increased efficiency is achieved by a 
dynamic programming approach exploiting the recursive 
structure of a closely related optimization problem. The 
bound for isotropic systems (Corollary [T]) is tight for the 
tested set of resources. It can be used to rule out distilla- 
tion of up to 9 fixed copies supporting the long standing 
conjecture of asymptotic impossibility of distillation in 
the isotropic case. The second bound (Corollary [5]) es- 
tablishes an efficient and generic connection between the 
distillable nonlocality of isotropic systems and the gen- 
eral case. It explains the increasing distillation potential 
when departing from the isotropic line, that has been 
observed in many distillation protocols, and is tight for 
correlated nonlocal boxes. 

The presented bounds are based on solutions to a re- 
cursive optimization problem that is new in this context. 
We believe that this abstraction can help to answer more 
general distillation questions and contributes to a deeper 
understanding of the possibilities of such protocols. It 
would be very interesting to find an efficient algorithm 
or an explicit solution or bound for the defined problem. 
This could be an important step towards proving general 
asymptotic (n —¥ oo) bounds to distillable nonlocality in 
the future. 
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Appendix A: A symmetry of distillation protocols 

Lemma 2. For any deterministic distillation protocol 
DP = (P n ,fo,fi,go,gi) H holds that 

NL (P n , f , h , 9o , 9i ) = NL (P n , /„ , h , g , g t ) 
= NL(P n J J 1 ,g ,g 1 ) 
= NL(P n J J 1 ,g ,g 1 ) 
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with f x {o) = 1 — fx{o) and g y (a) = 1 — g y (a) for all 
a G {0,1}" and all x, y G {0,1}. 

Proof. Observe that for any f,g G {0,1} 2 and any 
wiring distribution Was on n systems we have 

(f,g) = (l-2f) T W A B(l-2g) 

= -(1-2(1- f)) T W AB (l-2g) 
= -(l-2f) T W AB (l-2g) 
= -{f,9)- 

Apply it on 

NL(DP) = \(fo,g ) + (/o.ffi) + </i,<?o) - (h,9i)\ 
to prove the lemma. □ 



Appendix B: Calculate <5+ and S n by dynamic 
programming 

Lemma 3. Given numbers n > 1 and < k,l < 2 n , 

the values 5^(k,l) and S~(k,l) can be calculated in 2°^ 
steps. 

Proof. Taking advantage of the recursive structure of 5£ 
and 5~ we build the two tables 



D n = {$£{iJ)}o<i,j<2» and D n = {5 n (i,j)} < 



<i,j<2" 



recursively from the tables D^_ x and D n _ 1 . By read- 
ing out from D£_x and D~_ t the best of all possible 
solutions on for every entry D+(i,j) — 5^(i,j) 

and D~(i,j) = S~(i,j) we fill the tables £)+ and D~ 
element-wise. To save time some simple reduction rules 
derivable from the inner product definition and the triv- 
ial symmetries S+(i,j) = S+(j,i) and S~(i,j) = 6~(j,i) 
for each i, j can be used here - only an eighth of the 
table must actually be calculated. The base tables are 
given by DJ (k, I) = 1 and D^(k, I) = 1. Once the tables 
are filled, for fixed n and P the value D(n, P) can be 
bounded with the data in D+ and D~ by iterating over 
all possible Boolean function classes < |/ a T 1 (l)| < 2" 
and < l^^ 1 (1)| < 2". It requires roughly f™ =2 2 4l ~ 5 G 
2°C"0 steps to fill the tables in this way - a workload that 
is feasible for n < 10 on a normal PC. □ 



Appendix C: Finding the optimal Pi SO 

Bounding D(n,P) by Corollary [2] requires finding the 
least nonlocal isotropic system P iso such that P can be 
obtained from a convex combination of Pi SO and a local 
system Pi,. 



Lemma 4. Given a system P one can efficiently deter- 
mine the isotropic system Pi so , such that NL(Pi so ) is 
minimal under the constraint P = qPi so + Pl(1 — q) for 
any q G [0, 1] and any local system P F . 
Proof. The system P identifies a particular CHSH in- 
equality, now called CHSHp, and therefore a family of 
isotropic systems denoted 

P lso (e) = eP NL + (1 - e)P F 

with e G [0, 1] and Pnl as the nonlocal vertex corre- 
sponding to CHSHp and the unique isotropic system 
Pf G CHSHp. To find the optimal e we first calculate 
the local part of P [l4| with the linear program defined 
by Fitzi et al. (HI- The local part of P is the optimal 
value of the following linear program: 

max J2iPi 

SUCh that Y,iPi P L,iA V B( a ' b ) ^ P AB( a ' b ) 

Pi > 0. 

Here, Pl,% denote the vertices of the local polytope. Since 
Y^iPi i s maximal and there is only one nonlocal vertex 
violating CHSHp we have 



p=Y,PiPt,i + {i 



,Pi)p> 



NL- 



Now, we decompose the local system P* = J2j -Plj 
into 

P* =p f P F + (l- Pf )P L 

where Pl is any local system. By the local version of 
Lemma 1 in [15[ (holds since normalization and locality 
are linear properties) the maximal pj possible is given by 

a,b,x,y P F / B (a, 0) 

Therefore, the parameter e fixing the optimal Pi SO (e) can 
easily be constructed via 

p = Y,p 1 p* + (±-J2pd PNL 

i i 

= ^PiipfPp + (l-Pf)P L ) + (1 - £>)PjVL 

i i 

= ^Pi(l-p f )P L + Y,P*Pf P F + (1 - ^2pi)PNL 

i i i 

= ^2 Pi (l-p f )P L + (1 - 5> ( (1 -pf))P iS o(e), 

i i 

with 

1 -J2^P^ 



e = 



T,iPiPf + 1 -E l P» 

So, e can be determined efficiently. 



□ 
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